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Abstract 

Lepton flavor violation (LFV) is studied within a realistic unified framework, based on supersymmetric 
SO(10) or an effective G(224) = SU(2) L x SU(2) R x SU{4) C symmetry, that successfully describes (i) 
fermion masses and mixings, (ii) neutrino oscillations, as well as (iii) CP violation. LFV emerges as an 
important prediction of this framework, bringing no new parameters, barring the few SUSY parameters, 
which are assumed to be flavor-universal at M* > Mqut- We study LFV (i.e. /j, — > &y, r — ► /17, r — > &y 
and (J,N — > eiV) within this framework by including contributions both from the presence of the right 
handed neutrinos as well as those arising from renormalization group running in the post-GUT regime 
(M* — > Mgut)- Typically the latter, though commonly omitted in the literature, is found to dominate. 
Our predicted rates for /1 — > 67 show that while some choices of (m , mi/a) are clearly excluded by the 
current empirical limit, this decay should be seen with an improvement of the current sensitivity by a 
factor of 10-100, even if sleptons are moderately heavy (< 800 GeV, say). For the same reason, fj, — e 
conversion (fiN — > eTV) should show in the planned MECO experiment. Implications of WMAP and 
(g — 2) M -measurements are noted, as also the significance of the measurement of parity-odd asymmetry 
in the decay of polarized fi + into e +r y. 
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1 Introduction 



Individual lepton numbers (L e , L M and L T ) being symmetries of the standard model 
(SM)(with ml = 0), processes like /x — > e7, r — ► 7x7 and r — * e7 are forbidden within this 
model. Even within simple extensions of the SM (that permit m* 7^ 0), they are too strongly 
suppressed to be observable. Experimental searches have put upper limits on the branching 
ratios of these processes: Br(/x — > ej) < 1.2 x 10~ n [1], Br(r — > /X7) < 3.1 x 10~ 7 [2] and 
Br(r — > e7) < 3.7 x 10~ 7 [3]. The extreme smallness of these branching ratios poses a chal- 
lenge for physics beyond the standard model, especially for supersymmetric grand unified 
(SUSY GUT) models, as these generically possess new sources of lepton flavor violation that 
could easily lead to rates even surpassing the current limits. 

In this paper, we study how lepton flavor violation (LFV) gets linked with fermion masses, 
neutrino oscillations and CP violation within a predictive SUSY grand unified framework, 
based on either SO(10) [4], or an effective (presumably string derived) symmetry G(224) 
= SU(2) L x SU(2) R x SU(A) C [5]. The desirability of having an effective symmetry that 
possesses SU(4)-color [5], in view of the observed neutrino oscillations and the likely need 
for leptogenesis as a means for baryogenesis [6,7], has been stressed elsewhere (see e.g. [8]). 

A predictive framework based on supersymmetric SO(10) or G(224)-symmetry has been 
proposed by Babu, Pati and Wilczek (BPW) in [9], which successfully describes the masses 
and mixings of all fermions including neutrinos. In particular it makes seven predictions, all 
in good accord with observations. This framework was recently extended by us in Ref. [10] 
to describe the observed CP and flavor violations by allowing for phases in the fermion mass 
matrices. Remarkably enough, this extension could successfully describe the masses of all the 
quarks and leptons (especially of the two heavier families), the CKM elements, the observed 
CP and flavor violations in the K° — K° system (yielding correctly Am^ and ex) and the 
B^ — B^ system (yielding the correct values of Am^ and S^ s ). 

In this paper, we study lepton flavor violating processes, i.e. /x — > ej, r — * /17, r — > and 
fiN — > eN, within the same framework [9,10]. The subject of LFV has been discussed widely 
in the literature within supersymmetric extensions of the standard model. (For earlier works 
see Ref. [11-13]). Our work based on SUSY SO(10) or G(224) differs from those based on 
either MSSM with right-handed neutrinos (RHN's) [11,13,14] or SUSY SU(5) [15], because 
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for these latter cases the RHN's are singlets and thereby their Yukawa couplings are a 
priori arbitrary. By contrast, for G(224) or SO(10) the corresponding Yukawa couplings are 
determined in terms of those of the quarks at the GUT-scale (such as h(z/ T ) Dirac h top ) (see 
Ref. [9]). Thus the SUSY G(224)/SO(10)-framework is naturally more predictive than the 
MSSM or SUSY SU(5)-framework. 

In addition, our work differs from all others, including those based on SUSY SO(10) [16] 
as well, in two other important respects: First, we work within a predictive and realistic 
framework [9, 10] which (as mentioned above) successfully describes a set of phenomena - 
i.e. (a) fermion masses, (b) CKM mixings, (c) neutrino oscillations, (d) observed CP and 
flavor violations in the K and B systems, as well as (e) baryogenesis via leptogenesis [8]. As 
we will see, lepton flavor violation emerges as an important prediction of this framework, 
bringing no new parameters (barring the few flavor-universal SUSY-parameters). 

Second, we do a comprehensive study of LFV processes by including contributions from 
three different sources: (i) the sfermion mass-insertions, b % [ L RR , arising from renormalization 
group (RG) running from M* to M GUT ~2x 10 16 GeV (where M* denotes the presumed 
messenger-scale, with Mqut < M* < M string , at which flavor-universal soft SUSY breaking 
is transmitted to the squarks and sleptons, like in a mSUGRA model [17]), (ii) the mass- 
insertions (S l l L ) RHN arising from RG running from Mqut to the right handed neutrino mass 
scales M^, and (iii) the chirality-flipping mass-insertions 8 l [ RRL arising from A— terms that 
are induced solely through RG running from M* to M GUT involving gauginos in the loop. All 
the three types of mass-insertions: S l / L RR , (S l i L ) RHN and S % i RRL are in fact fully determined 
in our model. (See [10] for details). Most previous works in this regard have included 
only the second contribution associated with the RH neutrinos in their analysis. 1 We find, 
however, that it is the first and the third contributions associated with post-GUT physics 
that typically dominate over the second in a SUSY unified framework. 

A brief review of our previous work and our results are presented in the following sections. 

1 Barbieri, Hall and Strumia (in Ref. [12]) have discussed the relevance of the contributions from the mass-insertions 5^ L RR 
and those from the induced A— terms, but without a realistic framework for light fermion masses and neutrino oscillations. 
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2 A brief review of the BPW framework and its extension 



The Dirac mass matrices of the sectors u, d, I and v proposed in Ref. [9] in the context of 
SO (10) or G(224)-symmetry have the following structure: 
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a — e 





ff + e 
1 



M°- 



Ma 





r]' — e' 
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These matrices are defined in the gauge basis and are multiplied by ^ L on left and 
^/r on right. For instance, the row and column indices of M u are given by (uL,CL,ii) 
and (ur, cr, tfi) respectively. These matrices have a hierarchical structure which can be 
attributed to a presumed U(l)-flavor symmetry (see e.g. [8,10]), so that in magnitudes 1 
^> a ~ rj ~ e ^> Q2 ~ C22 ^ v' > e ' ■ The entries e and e' are proportional to B — L and are 
antisymmetric in family space (see below). Thus (e, e')— > — 3(e, e') as g — > /. Following the 
constraints of SO(10) and the U(l)-flavor symmetry, such a pattern of mass-matrices can be 
obtained using a minimal Higgs system consisting of 45h, 16h, 16h, 10h and a singlet S of 
SO(10), which lead to effective couplings of the form [8,10]: 

Auk = ^3316316310^ + [h 23 16 2 16 3 10 H (S/M) 

+ a 23 16 2 16 3 10 H (45 H /M r )(S/M) p + g 23 lQ 2 lQ 3 lQ d H (lQ H /M")(S/MY] 

+ [h 22 1 6 2 1 6 2 1 H (S/M) 2 + g 22 lQ 2 lQ 2 lQ d H (lQ H /M")(S/M) q+l ] 

+ [g l2 lQ l lQ 2 lQ d H {\6 H /M")(S/My +2 + a 12 16 1 16 2 10 H (45 H /M')(S/M) p+2 ] . (2) 

The mass scales M', M" and M are of order M str i n g or (possibly) of order Mqut [18]. 
Depending on whether M'(M") ~ Mgut or M str i n g (see [18]), the exponent p(q) is either one 
or zero [19]. The VEVs of (45 H ) (which is along B - L), (16 H ) = (16 h) (along (v RH )) and 
(S) are of the GUT-scale, while those of (10#) and (16#) are of the electroweak scale [9,20]. 
The combination 10 H -45 H effectively acts like a 120 which is antisymmetric in family space 
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and is along B—L. The hierarchical pattern is determined by the suppression of the couplings 
by appropriate powers of Mgut/{M, M'orM"). For details on how Eq. (J2J) emerges from 
Eq. (P) see Refs. [8-10]. 

The right-handed neutrino masses arise from the effective couplings of the form [21]: 

£ Maj = ^1646,16^16^ (3) 

where the fi/s include appropriate powers of (S)/M. The hierarchical form of the Majorana 
mass- matrix for the RH neutrinos is [9]: 



x z 
y 
z y 1 



M R (4) 



Following flavor charge assignments (see [8]), we have 1 ^> y ^> z ^> x. The magnitude 
of M R is estimated (with / 33 » 1, (16 H ) « 2 x 10 16 GeV and M » M st « 4 x 10 17 GeV) 
as [8,9]: M R = f 33 (16 H ) 2 /M « (10 15 GeV)(l/2-2). 

Thus the Majorana masses of the RH neutrinos are given by [8,9]: 

M 3 ^ M R ^ 10 15 GeV (1/2-1), M 2 « |y 2 |M 3 w (2.5xl0 12 GeV) (1/2-1), 

( 5 ) 

Mi w | a; - z 2 |M 3 ~ 10 10 GeV(l/4-2). 

Note that both the RH neutrinos and the light neutrinos have hierarchical masses. 

In the BPW model of Ref. [9], the parameters <J,t],e etc. were chosen to be real. To 
allow for CP violation, this framework was extended to include phases for the parameters 
in Ref. [10]. Remarkably enough, it was found that there exists a class of fits within the 
SO(10)/G(224) framework, which correctly describes not only (a) fermion masses, (b) CKM 
mixings and (c) neutrino oscillations [8,9], but also (d) the observed CP and flavor violations 
in the K and B systems (see Ref. [10] for the predictions in this regard). A representative of 
this class of fits (to be called fit A) is given by [10]: 

o- = 0.109-0.012i, rj = 0.122 - 0.0464i, e = -0.103, rj' = 2.4 x 10" 3 , 

(6) 

e' = 2.35 x 10"V( 69 °), C22 = 9.8 x lQ- 3 e -< 1490 \ {M° u , M° d ) « (100, 1.1) GeV. 
In this particular fit (% 2 is set to zero for the sake of economy in parameters. However, 
allowing for Q2 ~ (1/3) (C22) w °uld still yield the desired results. Because of the success 
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of this class of fits in describing correctly all four features (a), (b), (c) and (d)-which is a 
non-trivial feature by itself - we will use fit A as a representative to obtain the mass- insert ion 
parameters 8 l l L RR , (5^ L ) RHN and 8 l [ RRL in the lepton sector and thereby the predictions of 
our model for lepton flavor violation. 

The fermion mass matrices M u , Md and Mi are diagonalized at the GUT scale »2x 10 16 
GeV by bi-unitary transformations: 

KZ = X L ,m M uAl X^ Al) (7) 

The analytic expressions for the matrices Xf jR can be found in [10]. The corresponding 
expressions for X l L R can be obtained by letting (e, e')— > — 3(e, e'). 
We now discuss the sources of lepton flavor violation in our model. 



3 Lepton Flavor Violation in the SO(10)/G(224) Framework 

We assume that flavor-universal soft SUSY-breaking is transmitted to the SM-sector at a 
messenger scale M*, where M GUT < M* < M string . This may naturally be realized e.g. 
in models of mSUGRA [17], or gaugino-mediation [22]. With the assumption of extreme 
universality as in CMSSM, supersymmetry introduces five parameters at the scale M*: 

m o, ™>i/2, A>, tan/3 and sgn(/i). 

For most purposes, we will adopt this restricted version of SUSY breaking with the added 
restriction that A Q = at M* [22] . However, we will not insist on strict Higgs-squark-slepton 
mass universality. Even though we have flavor preservation at M*, flavor violating scalar 
(mass) ^transitions arise in the model through RG running from M* to the EW scale. As 
described below, we thereby have three sources of lepton flavor violation. 
(1) RG Running of Scalar Masses from M* to M GUT . 

With family universality at the scale M*, all sleptons have the mass m c at this scale and 
the scalar (mass) 2 matrices are diagonal. Due to flavor dependent Yukawa couplings, with 
h t = h b = h T (= /i 33 ) being the largest, RG running from M* to M G tjt renders the third 
family lighter than the first two (see e.g. [12]) by the amount: 

SOm 2 

Amf L = A ™t = A < = = A w " ( 16^) h " HM*/M GUT ) • (8) 
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The factor 30^12 for the case of G(224). The slepton (mass) 2 matrix thus has the form 
= diag(m 2 , m 2 , m 2 — A). As mentioned earlier, the spin-1/2 lepton mass matrix is 
diagonalized at the GUT scale by the matrices X l L R . Applying the same transformation 
to the slepton (mass) 2 matrix (which is defined in the gauge basis), i.e. by evaluating 
X^(M~°^)ii Xl and similarly for L^R, the transformed slepton (mass) 2 matrix is no longer 
diagonal. The presence of these off-diagonal elements (at the GUT-scale) given by: 



(9) 



l 3 



induces flavor violating transitions l l L R — > V L R . Here m z ~ denotes an average slepton mass 
and the hat signifies GUT-scale values. 

(2) RG Running of the A— parameters from M* to Mqut- 

Even if A D = at the scale M* (as we assume for concreteness, see also [22]). RG running 
from M* to Mqut induces A— parameters at Mqut, invoving the SO(10)/G(224) gauginos; 
these yield chirality flipping transitions {l l L)R — > l RL )- 

Evaluated at the GUT scale, the A— parameters, induced respectively through the cou- 
plings hij, ciij and gij, are given by: 

4 } = f ^glMMiE-) (iJ = 2, 3) 



A 



(2) 



95 1 Oi J -(45 H > 2 Juf 
2 8tt 2 M' yiO iW A m VMi (H 

9ij{Wn) n 2 M ^ 7rl/ 

2 8tt 2 



{ij = 23, 12) 



(10) 



4? = f^^gfoMMtE-) (ij = 23, 22, 12) 



The coefficients ^, ^) are the sums of the Casimirs of the SO(10) represen- 
tations of the chiral superfields involved in the diagrams. For the case of G(224), 
(f , f , f )->(§ , f , f ). Thus, summing A<U, A^ and A (3) , the induced A matrix for the 



leptons is given by: 
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Z t K 
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(11) 



where Z = (^)h t gf M x , K w = ln(j^) and K 16 = In(j^). For simplicity if we let 
Mi 6h ~ Mio H ~ M G ut, we can write the A matrix in the SUSY basis as: 

-285e' + 90?/ 



A 1 

Si 1 



LR 



M* 



UT 



285e' + 90r/' 90C| 2 - 27C 2 M 2 -285e + 90/7 - 27a 
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X l R (l2) 
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Approximate analytic forms for R are given in Ref. [10], and X^ R can be obtained from 
X-l,r by the substitutions (e, e') — > — 3(e, e'). The chirality flipping transition angles are 
defined as : 

= ( Al LRh (^) = {A R )ij ( tan/3 U m 2) • ( 13 ) 

(3) RG Running of scalar masses from M G ut to the RH neutrino mass scales: 

We work in a basis in which the charged lepton Yukawa matrix and are diagonal 
at the GUT scale. The off-diagonal elements in the Dirac neutrino mass matrix Y^r in this 
basis give rise to lepton flavor violating off-diagonal components in the left handed slepton 
mass matrix through the RG running of the scalar masses from Mqut to the RH neutrino 
mass scales Mj^. The RH neutrinos decouple below Mr.. (For RGEs for MSSM with RH 
neutrinos see e.g. Refs. [13] and [23].) In the leading log approximation, the off-diagonal 
elements in the left-handed slepton (mass) 2 -matrix, thus arising, are given by: 

(OS™ = ~ {Sm lt Al) M^f 1 ) ■ (i4) 

71 k=l Rk 

The superscript RHN denotes the contribution due to the presence of the RH neutrinos. 
We remind the reader that the masses M^. of RH neutrinos are well determined within our 
framework to within factors of 2 to 4 (see Eq. (J3J). The total LL contribution is thus: 

W = + (S 1 llW N (15) 

Now, most authors including those using SUSY SU(5) with RHN's or SUSY SO(10) 
[15, 16] have cosidered only the second term (S l LL ) RH that arises due to the right-handed 
neutrinos. As mentioned in the introduction, however, the first term S l LL and the contribution 
of the A— term S l LRRL (Eq. (JEH)) are found to dominate over the second term (as long as 
\ii(M*/Mgut) ~ 1)- We obtain our results by including the contributions from all three 
sources listed above in Eqs. Q, (HHJ) and (fTlj) . They are presented in the following section. 

4 Results 

The decay rates for the lepton flavor violating processes U — > L/7 {i > j) are given by: 

2 ^ 

ra + -/;7) = ^(l^| 2 + K| 2 ) (16) 
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Here A°l is the amplitude for — > (lj) +/ y decay, while A 3 R = A((l i )^ i — > (lj) +/ y). 
The amplitudes A 3 £ R are evaluated in the mass insertion approximation using the 
(^ll) T °*) 3rri ^lrrl calculated as above. The general expressions for the amplitudes A°l R 
in one loop can be found in e.g. Refs. [13] and [23]. We include the contributions from both 
chargino and neutralino loops with or without the /z— term. 

We evaluate the amplitudes by first going to a basis in which the chargino and the 
neutralino mass matrices are diagonal. Analytic expressions for this diagonalization can be 
obtained in the approximation \M% ± fi\ and |§Mi ± /i| ^> rriz and |M2/i| > m 2 ^ sin 2/3 [24]. 
This approximation holds for all the input values of (m , mi/ 2 ) that we consider. 

In Table 1 as well as in Fig. 1, we give the branching ratios of the processes /i — > 
ej, r — > fij, t — > e7 for the case of SO(10), with some sample choices of (m D , m 1( / 2 ). 
For these calculations we set ln^-^g^) = 1, i.e. M* « 3M GUT , tan/3 = 10, M Rl = 10 10 
GeV, M R2 = 10 12 GeV and = 5 x 10 14 GeV (see Eq. ©), and A a ( at M*) = 0. The 
corresponding values for G(224) are smaller approximately by a factor of 4 to 6 in the rate, 
provided ln(M* / 'Mqut) is the same in both cases (see comments below Eqs. (jS)) and ()10j0. 



(m , mi/ 2 )//tan/3 


Br(/i — > ej) 


Br(r -> 


Br(r — > ej) 




n > 


/i < 


n > 


/i < 


/i > 


/i < 


I (600, 300)//10 


3.3xl0~ 12 


9.8xl0~ 12 


2.4xl0~ 9 


3.1 xl0~ 9 


2.4xl0~ 12 


3.3xl0~ 12 


II (800, 250)//10 


2.9xl0~ 13 


1.7xl0~ 12 


1.9xl0~ 9 


1.9xl0~ 9 


2.0xl0~ 12 


2.0xl0~ 12 


III (450, 300)//10 


2.7xl0- n 


4.6xl0- n 


2.7xl0~ 9 


5.6xl0~ 9 


2.7xl0~ 12 


6.1xl0~ 12 


IV (500, 250)//10 


5.9xl0^ 12 


1.9xl0^ n 


4.8xl0~ 9 


6.4xl0" 9 


5.0xl0~ 12 


6.9xl0~ 12 


V (100, 440)//10 


1.02xl0~ 8 


1.02xl0~ 8 


8.3xl0~ 8 


8.4xl0~ 8 


1.0xl0~ 10 


1.0xl0~ 10 


VI (1000, 250)//10 


1.6xl0~ 13 


5.6xl0~ 12 


9.5xl0~ 10 


9.0xl0~ 10 


l.OxlO^ 12 


9.5xl0~ 13 


VII (400, 300)//20 


9.5xl0" 12 


3.8xl0- n 


1.4xl0~ 8 


1.8xl0" 8 


1.5xl0- n 


1.9xl0- n 



Table 1. Branching ratios of U — * Ijj for the SO(10) framework with k = ln(M* / 'Mgut) = 1; (m , m^) 
are given in GeV, which determine /1 through radiative electroweak symmetry breaking conditions. The entries for 
Br(/i — > ey) for the case of G(224) would be reduced by a factor f»4-6 compared to that of SO(10) (see text). 



To give the reader an idea of the magnitudes of the various contributions, we exhibit in 
table 2 the amplitudes for the process /i — > calculated individually from the four sources 
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8% L , 5 j i R)RL and (5 j ^ L ) RHN (see Eqs. ©, (JH and d), for a few cases of table 1. 



(m , m 1/2 )(GeV) 




A { ?(6lr) 


A R (S RL ) 




I, (600, 300) 


3.3 x 10~ 13 


-6.7 x 10~ 13 


-5.9 x 10- 13 


2.4 x 10" 14 


II, (800, 250) 


2.9 x 10~ 13 


-1.8 x 10- 13 


-1.6 x 10~ 13 


2.0 x 10" 14 


IV, (500, 250) 


4.8 x 10~ 13 


-9.7 x 10~ 13 


-8.5 x 10~ 13 


3.4 x 10~ 14 



Table 2. Comparison of the various contributions to the amplitude for fi —* e'y for cases I, II and IV, with /i > 0. 
Each entry should be multiplied by a common factor a . Imaginary parts being small are not shown. Note that 
columns 2,3 and 4 arising from RG running from M* — > Mgut (see text) dominate over the RHN contribution. 

Glancing at tables 1 and 2, the following features of our results are worth noting: 

(1) It is apparent from table 2 that the contribution due to the presence of the RH neutrinos 2 
(fifth column) is about an order of magnitude smaller, in the amplitude, than those of the 
others (proportional to 5 % [ L , 5 % l R and S RL ), listed in columns 2, 3 and 4. The latter arise 
from RG running of the scalar masses and the A— parameters in the context of SO(10) or 
G(224) from M* to Mgut- It seems to us that the latter, which have commonly been omitted 
in the literature, should exist in any SUSY GUT model for which the messenger scale for 
SUSY-breaking is high (M* > Mgut), as in a mSUGRA model. The inclusion of these new 
contributions to LFV processes arising from post- GUT physics, that too in the context of a 
predictive and realistic framework, is the distinguishing feature of the present work. 

(2) Again from table 2 we see that the two dominant contributions to At, = A(fi^ — > e + 7), 
arising from 8ll and c^^-insertions, partially cancel each other if /i > 0; they would however 
add if n < 0. By contrast, A R gets contribution dominantly only from 8 R l (column 4). 3 As 
a result we find that in our model, typically, \A R \ > \Al\ if /x > and \Al\ > \Ar\ if // < 0. 

(3) Owing to the general prominence of the new contributions from post-GUT physics, we 

2 In the context of contributions due to the RH neutrinos alone, there exists an important distinction (partially observed by 
Barr, see Ref. [16]) between the hierarchical BPW form [9] and the lop-sided Albright-Barr (AB) form [25] of the mass-matrices. 
The amplitude for fi — > cy from this source turns out to be proportional to the difference between the (23)-elements of the 
Dirac mass-matrices of the charged leptons and the neutrinos, with (33)-element being 1. This difference is (see Eq. Q) is 
r] — a fs; 0.041, which is naturally small for the hierarchical BPW model (incidentally it is also V c j,), while it is order one for 
the lop-sided AB model. This means that the rate for fj, —* ej due to RH neutrinos would be about 600 times larger in the AB 
model than the BPW model (for the same input SUSY parameters). 

3 Although &b,r 'S comparable to 8ll, its contribution to Ar (via the bino loop) is typically suppressed compared to that of 
S^i, to Al (in part by the factor (ol\ / 'a2)(Mi / 'M2)) in most of the parameter space. 
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see from table 1 that case V, (with low m and high m^) is clearly excluded by the empirical 
limit on /x — > e7-rate (see Sec. 1). Case III is also excluded, for the case of SO(10), yielding 
a rate that exceeds the limit by a factor of about 2 (for k = ln(M* / M GUT ) > 1), though 
we note that for the case of G(224), Case III is still perfectly compatible with the observed 
limit (see remark below table 1). All the other cases (I, II, IV, VI, and VII), with medium 
heavy (~ 500 GeV) to moderately heavy sleptons (800-1000 GeV), are compatible with the 
empirical limit, even for the case of SO(10). The interesting point about these predictions 
of our model, however, is that /i — > &y should be discovered, even with moderately heavy 
sleptons, both for SO(10) and G(224), with improvement in the current limit by a factor of 
10-100. Such an improvement is being planned at the forthcoming MEG experiment at PSI. 

(4) We see from table 1 that r — > /ry (leaving aside case V, which is excluded by the limit 
on fi — ► ej), is expected to have a branching ratio in the range of 2 x 10~ 8 (Case VII) to 
about (1 or 2) x 10~ 9 (Case VI or II). The former may be probed at BABAR and BELLE, 
while the latter can be reached at the LHC or a super B factory. The process r — > cy would, 
however, be inaccessible in the foreseeable future (in the context of our model). 

(5) The WMAP-Constraint: Of the cases exhibited in table 1, Case V (m = 100 
GeV, mi/2 = 440 GeV) would be compatible with the WMAP-constraint on relic dark 
matter density, in the context of CMSSM, assuming that the lightest neutralino is the LSP 
and represents cold dark matter (CDM), accompanying co-annihilation mechanism. (See 
e.g. [26]). As mentioned above (see table 1), a spectrum like Case V, with low m and higher 
mi/2, is however excluded in our model by the empirical limit on /i — > ej. Thus we infer that 
in the context of our model CDM cannot be associated with the co- annihilation mechanism. 

Several authors (see e.g. Refs. [27] and [28]), have, however considered the possibility 
that Higgs-squark-slepton mass universality need not hold even if family universality does. 
In the context of such non-universal Higgs mass (NUHM) models, the authors of Ref. [28] 
show that agreement with the WMAP data can be obtained over a wide range of mSUGRA 
parameters. In particular, such agreement is obtained for (m^,/m ) of order unity (with 
either sign) for almost all the cases (I, II, III, IV, VI and VII) 4 , with the LSP (neutralino) 
representing CDM. 5 (Here = sign{m 2 Hu d )\J\ m2 H u J> see [28]). All these cases (including 

4 We thank A. Mustafayev and H. Baer for private communications in this regard. 

5 We mention in passing that there may also be other posibilities for the CDM if we allow for either non-universal gaugino 
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Case III for G(224)) are of course compatible with the limit on fi — > cy. 

(6) Coherent /i — e conversion in nuclei: In our framework, fi — e conversion (i.e. 
+ N — > e _ + iV) will occur when the photon emitted in the virtual decay /i — > e7* is 

absorbed by the nucleus (see e.g. [29]). In such situations, there is a rather simple relation 
connecting the /i — e conversion rate with B(/i — > cy): B(fi — > ej) / (ou conversion / u ca pture) = 
R ~ (230 — 400), depending on the nucleus. For example, R has been calculated to be 
R ~ 389 for 27 Al, 238 for 48 Ti and 342 for 208 Pb in this type of models. (These numbers 
were computed in [29] for the specific model of [12], but they should approximately hold 
for our model as well.) With the branching ratios listed in Table 1 (~ 10 -11 to 10~ 13 ) for 
our model, oo convers i on / uu capt ure — (40-1) xl0~ 15 . The MECO experiment at Brookhaven is 
expected to have a sensitivity of 10~ 16 for this process, and thus will test our model. 

(7) Parity odd asymmetry in /x+ — > e+7 decay: Parity violation can be observed 
by studying the correlation between the momentum p e of e+ in ji + — > e+7 decay and 

— * 

the polarization vector P of positive muons (from tc + decays). The distribution of e+ 
is proportional to (1 + A p e .P) where A is the P-odd asymmetry parameter given by 
A((i + -> e+7) = (\A L \ 2 - \A R \ 2 )/(\A L \ 2 + \A R \ 2 ). Here A L is the amplitude for fj+ -> e+7 
decay, while Ar = A{^ — > e+7). In our model, as noted in (2), we typically have 
\A R \ > \A L \ and thus A(fx + — >■ e+7) < if fx > 0, and \A L \ > \A R \ and thus A > if 
pi < 0. For example, with (ra , mi/ 2 ) = (800, 250) GeV, fx > and tan/3 = 10, we obtain 
|A L | = |^ 1} (5 LL ) + Af(5 Li? ) + Af| = 1.3 x 10- 13 (see table 2) while ~ 1.6 x lO" 13 , 
and thus A ~ —0.25, while for (m ,m 1 / 2 ) = (500, 250) GeV and tan/3 = 10 we get, 
\A L \ ~ 4.7 x 10~ 13 and ~ 8.6 x 10~ 13 , yielding A ~ —0.54. The precise prediction of 
our model for A would thus be definitive once the SUSY spectrum is known. 

We can compare the predictions of our model for A with those of other SUSY models. 
In the MSSM with v R , since LFV arises through S LL type mixings, A L 3> A R , and thus 
A(fi + — > e+7) « +1, at least for tan/3 < 30 or so, regardless of the choice of (m , TOi/ 2 )- In 
SUSY SU(5) GUT, with or without z/r, the GUT threshold effects realized in the regime 
Mqut < A 4 — ^* generate type mixings, and will lead to ^> and thus A ~ — 1. In 
the SUSY 5*0(10) models with symmetric mass matrices, such as the ones studied in [12,30], 
Al = A R from GUT threshold effects, leading to a vanishing A. Thus, we see that a 

masses, or axino or gravitino as the LSP, or R-parity violation (with e.g. axion as the CDM). 
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determination of A may help sort out the specific type of GUT that is responsible for LFV. 
(8) Correlation between muon g — 2 and /i — > cy: Currently there exists a discrepancy 
between theory and experiment in the anomalous magnetic moment of the muon: Aa M = 
a^ xpt — a^ M = 251(93) x 10~ n [31]. This is a 2.7 sigma effect 6 and may be an indication of 
low energy supersymmetry. In our framework, this discrepancy can be considerably reduced 
for some, but not all, choices of the SUSY spectrum. When the sleptons are relatively 
light (< 500 GeV) with tan (3 = 10 — 20, the SUSY contribution to a M is in the range 
(50 — 200) x 10~ n . For example, following a recent numerical analysis (see [32] and references 
there in), we find Aa^ USY rs 180 x 10~ n for the cases of both IV and VII (see table 1). Note 
that when the SUSY contributions to Aa M becomes significant, B(fi — > cy) is enhanced. 
Thus, a confirmation of new physics contribution to a M , for example by improved precision 
in the e + e~ — > hadron data and in the theoretical analysis, would imply (in the context of 
a SUSY-explanation) that /i — > is just around the corner, within our framework. 

In summary, lepton flavor violation is studied here within a predictive SO(10)/G(224)- 
framework, possessing supersymmetry, that was proposed in Refs. [9,10]. The framework 
seems most realistic in that it successfully describes five phenomena: (i) fermion masses 
and mixings, (ii) neutrino oscillations, (iii) CP violation, (iv) quark flavor-violations, as well 
as (v) baryogenesis via leptogenesis [7]. LFV emerges as an important prediction of this 
framework bringing no new parameters, barring the few flavor-preserving SUSY parameters. 

As mentioned before, the inclusion of contributions to LFV arising both from the pres- 
ence of the RH neutrinos as well as those from the post-GUT regime, that too within a 
realistic framework, is the distinguishing feature of the present work. Typically, the latter 
contribution, which is commonly omitted in the literature, is found to dominate. Our results 
show that - (i) The decay fj, — > should be seen with improvement in the current limit by a 
factor of 10 - 100, even if sleptons are moderately heavy (~ 800 GeV, say); (ii) for the same 
reason, /j, — e conversion (/iN — > eN) should show in the planned MECO experiment, and 
(iii) r — > /x7 may be accessible at the LHC and a super B-factory. It is noted that the muon 
(g — 2)-anomaly, if confirmed, would strongly suggest, within our model, that the discovery 
of the n — > e7 decay is imminent. The significance of a measurement of the parity-odd 

6 This analysis is based on theory and data on e+e~ — » hadron. If t — * v T + hadron data is used, this discrepancy reduces 
to 1.3 sigma; this may however be less reliable [31]. 
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asymmetry in polarized fi + decay into e + 7 is also noted. In conclusion, the SO(10)/G(224) 
framework pursued here seems most successful on several fronts; it can surely meet further 
stringent tests through a search for lepton flavor violation. 
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Fig. 1. Log of Br(/i — > e-y) divided by the experimental bound (1.2 x 10 11 ) obtained for the SO(10) framework 
with ln(M* /Mgut) = 1, tan/3 = 10 and fi > vs m D (in GeV) with m 1/2 = 200, 250 and 300 GeV. 
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